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Abstract. The subdwarf B star PG 1605+072, with an unusually low log (g/cms~ 2 ) 
-5.3, shows a high number of non-radial pulsation modes, making it a promising candi- 
date for asteroseismology. This could allow probing the star's interior to gain important 
p/ insights in its structure and evolution. Comparison of previous work conducted over 

rys. . the last decade shows clear amplitude variation and hints of frequency variation. We 

analyse white light photometric data of the Multi-Site Spectroscopic Telescope (MSST) 
i-C \ and Whole Earth Telescope (WET) XCov22 campaigns using prewhitening techniques 

♦~H- and O-C diagrams. A total of 85 significant frequencies are identified, among them 

q ' more than 20 frequency sums and harmonics. Moreover, it is shown that the main 

5h , mode's amplitude varies like a sine with a period of ~ 630 d, indicative of long-term 

t/5 ■ beating. Strong hints for the existence of frequency multiplets support the hypothesis 

& ' that PG 1605+072 is a slow rotator with V eq < 0.9 kms , contrary to previous claims. 

Existing asteroseismic models mainly suggest that the star possesses a high mass of 

~ 0.7 M Q , presuming the main pulsation mode to be radial. This calls for an alternate 

^ formation channel. 

OO 

t> 

^H 

ry-) ' 1. Introduction 

o" 

The subdwarf B star PG 1605+072 , also known as V338 Ser, was found to be a p- 
mode pulsator by iKoen et al.l (1 1998b . showing more than 30 pulsation frequencies. The 



main mode at 2076 mHz (~ 179 cd ) initially proved to be very strong, possessing 
an amplitude of more than 60 mmi. This rich pulsation spectrum wa s confirmed by 
the tw o-week photometric multi-site campaign that was conducted by Kilkenny et al.1 
(1 19991) in 1997, revealing a total of 55 pulsation frequencies, with periods ranging from 
200 s to 600 s. However, the authors found that the main mode's amplitude decreased 
to less than 30 mmi on t he timescale of roughly one year. The analysis of spectroscopic 
and photometric data by lO' Toole et al.l (12000). obtained in 1999, showed that the main 
mode not only weakene d, but completel y dis appeared. In spectr oscopic campaigns con- 
ducted in 2000 bv lQ' Toole etaD (120021) and lWoolf et"aD (120021) . the mode returned but 
wa s still weaker than in 1997. Finally, spectroscopic and photometric data from 2001 
by Falter et all (120031) produced a spectrum similar to the findings in iKilkenny et al. 



1999). This behaviour hints at amplitude variation on the timescale of years. 



Spectroscopically, PG 1605+072 shows T eS = 32 100 ± 1 000 K, which is typical 



of a p -mode pulsating sdB star, but also a low log (g/cm s~ 2 ) = 5.25 + 0.10 (IKoen et al. 
1 19981. confirmed by later observations), putting the star in an unusual position in the 
Hertzsprung-Russell diagram, about 0.6 dex away from the zero-age extreme horizontal 
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branch. This can either be explaine d by an evolved st atus (IDorman et al.lll993l) or an 
unusually high mass (see, e. g. Jvan Grootel et aLll2010l) . 

Anothe r prob lem which has not been solved to this day is the rotation velocity. 
iHeber et all d 19991) considered the line broadening observed in PG 1605+072 to be 
caused by rotation, delivering a rotation velocity of 39kms _1 . This value is unusu- 
ally high for an sdB star, especially considering that PG 1605+072 seems to be a single 
star: It will evolve to a single white dwarf and keep its high rotation velocity (if there 
is no mechanism to get rid o f the angular mom entum), but single white dwarfs are 
measu red to be slow rotators (IHeber et al.lll997l and references therein). iGeier et aL 
(1201 Ol) also measure rather low rotational velocities for sdB stars, namely ~ 8 km s 



for the majority of their targets. On the other hand, a high rotation velocity would be 
able to explain the rich pulsation spectrum, as it lifts the (21 + l)-fol d degene r acy of 
pulsation modes with the same k and / but different m values. Thus, iKawalerl (1 19991) 
interpreted the spectrum as mainly being the result of trapped non-radial modes. The 
authors estimated a rotation velocity of ~ 130 km s' 1 . 

However the line broadening can also be interpreted as pulsational broadening. 



Kuassivi et al.1 J2005I) came to the conclusion that the rapid Doppler shift of the lines 



due to pulsation got smeared in the data of IHeber et all because of the long integra- 
tion time of 600 s. Using FUSE (Far Ultraviolet Spectroscopic Explorer) data, they 
found a Doppler shift of ~ 17 km s -1 . The order of m agnitude is consistent with 
the results of b'Toole et all (l2002h and IWoolf et all (|2002|) as well as the results from 



the spectroscopic part of the Multi-Site Spectroscopic Telescope (MSST) campaign 
(lO'Toole et al.ll2005r) . The complexity of the pulsation spectrum could be explained 



by assuming the low-amplitude modes to be second- and third-order harmonics as well 
as non-linear combinati ons of the very strong main modes. Following this idea, the 
pulsation spectrum from Kilkenny et all (1 19991) can be built up using just 22 basic fre- 
quencies (Ivan Grootel et al.ll2010h . 

We investigate these problems using photometric data from the MSST and Whole 
Earth Telescope (WET) XCov22 campaigns (see Sect. |2]). In Sect. [3j the results from 
the prewhitening process are presented. Evidence for amplitude variation is followed in 
Sect.@]and mode splitting is analysed in Sect. [5] The implications for our understanding 
of PG 1605+072 are then discussed in Sect. [6] before the conclusions are drawn in 
Sect |7] 



2. Data 



PG 1605+072's rich frequency spectrum, high amplitudes and long pulsation periods 
were considered excellent prerequisites for detecting spectroscopic variations and ap- 
plying asteroseismic methods to understand its evolutionary status. Thus, in 2002, the 
Multi- Site Spectroscopic Telescope (MSST) campaign was conducted by IHeber et aL 
(2003) in order to obtain simultaneous high quality spectroscopic and photometric data, 
allowing the resolution of more pulsation frequencies without strong aliasing effects. 
Additionally, the MSST project was supported by the selection of PG 1605+072 as an 
alternate target in the Whole Earth Telescope (WET) XCov22 photometric campaign 

(ISchuh et alJl2003l) . 

Parts of the campaign data have already been analysed (lO'Toole et al.l l2004). We 
analysed the combined photometric data, spanning a time base of approximately 30 
days. The time stamps have been converted into BJDtdb- Since many individual light 
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curves from different sites overlapped in time, cross-correlation could be used to un- 
cover and if necessary eliminate remaining timing errors. Overlapping light curve frag- 
ments with timestamps known to be already in the correct format were used to create 
a model light curve in order to cover the whole time base. The remaining data could 
be compared to this model using cross-correlation again. This way, timing problems 
could be identified and treated individually for each light curve. After this procedure, a 
timing error larger than ~ 0.33 s is unlikely to persist. 

For the further analysis, two versions of the total light curve were used. First, the 
original, non-weighted light curve was used for prewhitening (see Sect. 13.11) . Then, 
weights were estimated calculating a moving standard deviation of the residual light 
curve. In the overlapping regions, only the light curve data set with the highest weights 
was retained. 



3. Prewhitening 

3.1. Non- weigh ted light curve 

Using the program Period©4 by P. Lenz and M. Breger, the original, non-weighted 
light curve was prewhitened until residual peaks fell below a limit of four times the 
noise level. First, the total light curve was divided into two roughly equal parts with 
a time basis of 16.0 d and 14.6 d respectively, leading to a frequency resolution A/ ~ 
0.062 cd _1 and ~ 0.068 cd _1 . These light curves were prewhitened individually in 
order to look for frequencies present in both halves of the data set. Due to the inho- 
mogeneity of the photometry, 75 frequencies could be fitted for the first half's data, but 



only 41 for the second half (see Langfellner et al., in preparation, or Langfellner l2011 
for details). As a result, only 30 frequencies were found to be present in both halves. 

Prewhitening of the total light curve with a frequency resolution of ~ 0.033 cd _1 
yielded 85 frequencies in total, situated in the range between 150 and 550cd _1 . For 
higher frequencies, no significant frequencies were spotted. For frequencies below 
~100cd -1 , pulsation frequencies become undistinguishable from atmospheric effects. 

Snapshots of the prewhitening process are shown in fig. \T\ The 10 frequencies 
with the highest amplitudes are listed in table [Q Very close to the main mode, another 
strong peak was detected. As the frequency difference is far smaller than the resolution 
and the second peak had not been detected in the split light curves, we assume this peak 
to be artificial and give it the rank "la". The combination of f\ and f\ a should result in 
an amplitude near 50mmi, as visible in fig. CD 

Out of the 85 frequencies, 18 frequency sums and the first harmonics of the 
strongest three modes could be found with deviations Sf well below the frequency res- 
olution. Remarkably, all combinations of the main mode with the other strong modes 
up to rank 16 are present in the data, apart from ranks 13 and 15. Furthermore, the 
three highest frequencies consist of the sums of the first harmonics of f\ or fi and an 



other strong frequency. This is indicative of non-linear pulsations (lAntonello & Aikawa 



1998). 



3.2. Comparison with Kilkenny et al. (1999) 

Let us compare the results with those from lKilkenny et all (11999). representing the most 



extensive previous photometric multi-site campaign on PG 1605+072. Fortunately, at 
least 32 of the total 55 frequencies match the MSST/WET results, i. e. they differ by 
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Figure 1. Prewhitening process of the original, non-weighted light curve. The 
uppermost panel shows the original Fourier spectrum, the panels below the residual 
spectra after prewhitening 10, 30 and 85 frequencies. 



less than A/. At least one additional frequency matc h could be ex plained by aliasing 



(lSO^ecd" 1 in MSST/WET and 181.475 cd _1 in iKilkennv et alJ data). However, 
there are also stronger frequencies in the MSST/WET data, such as f$, that did not 
have any corresponding frequency detection. 

As expect ed, the MS ST and WET campaigns, covering twice the time of the 
Kilkenny et alJ campaign, deliver more significant frequencies - 85 compared to 55. 
However, the number of frequencies confirmed by appearing in both halves of the data 
is smaller (22 vs. 28) due to the inhomogeneity of the MSST/WET data. Besides, the 
number of frequencies that can be constructed by the sum of two other frequencies has 
almost doubled - from 10 to 18. Additionally, the first harmonics can be found for the 
three strongest frequencies instead of just for the main mode. 

It is remarkable as well that the frequency range could be extended by identifying 
three significant frequencies in the region between 500 and 600 cd" 1 . Thus, the pulsa- 
tion periods re ach down at least to 160 s. However, higher harmonics had already been 
found by lKoen et alJ dl99a) . reaching up to the third harmonic of fi. 

Apart from the frequency values, the comparison of the amplitudes gives interest- 
ing results. The main mode is much stronger in the MSST/WET data, even if the two 
frequencies ranked "1" and "la" are thought to interfere destructively. This frequency 
splitting in the Fourier spectrum is another indication for amplitude variation. 
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Table 1. The ten strongest frequencies in MSST/WET non-weighted light curve. 
Note the splitting of the main mode f\ into two frequencies with ranks "1" and 
"la". The rank of /io is set in brackets, as this frequency could not be resolved in 
both halves of the light curve. See Langfellner et al. (in preparation) for a complete 
frequency table. 



Rank 


/ 


A 


P 


Freq. sum 


Sf 




[cd^ 1 ] 


[mi] 


[s] 




[cd 1 ] 


4 


163.415 


0.0101 


528.71 






3 


171.574 


0.0121 


503.57 






la 


179.343 


0.0205 


481.76 






1 


179.347 


0.0668 


481.75 






6 


181.608 


0.0067 


475.75 






2 


181.659 


0.0151 


475.62 






(10) 


181.703 


0.0055 


475.50 






8 


236.997 


0.0062 


364.56 






5 


238.752 


0.0081 


361.88 






9 


245.781 


0.0058 


351.53 






7 


350.920 


0.0062 


246.21 


/i+/3 


0.001 



Another remarkable result is the freque ncy trip let comprising fa, f% and /io (see 
tabled]), as it is also identifiable in the lKilkenny et all data, but with completely different 
amplitudes. In general, the amplitude values change a lot - a frequency ranked as high 
as "5" in the MSST/WET data even does not appear at all in the old data. 

3.3. Using weights 

As mentioned in Sect. El the residuals of the non-weighted prewhitening process were 
used to derive point weights that served as a criterion for choosing the best light curve 
fragment in overlapping regions. For the weight derivation, a Fortran program by 
Roberto Silvotti was used that calculates a moving standard deviation of the residuals. 

This modified, weighted light curve was used to repeat the prewhitening pro- 
cess. As before, 85 frequencies could be identified. Taking the frequency resolution 
into account, we found that 71 of the frequencies are identical to the ones from the 
non-weighted prewhitening procedure. Four more frequency sums could be identified, 
yielding a total of 22. Again, we found the first harmonics of the three strongest modes. 

Considering the main mode, its double feature disappeared and showed only a 
single peak with an amplitude reduced to ~ 51 mmi. The amplitudes of the other 15 
strongest frequencies only changed slightly (i. e. < 10%). Merely for amplitude ranks 
> 20, particularly > 60, frequencies that were not found in the non-weighted light 
curve were spotted. In general, however, these differences are not substantial. Conse- 
quently, the non-weighted light curve could be used for the further analysis. 



4. Amplitude variation 

For the search for amplitude variation visible on the timescale of the MSST/WET obser- 
vations, the O-C diagram with its high sensitivity proved to be the tool of choice. For 
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Figure 2. Long-term beating of main mode f\. Together with the 40 amplitude 
values from O-C diagrams for f\ in the MSST/ WET data (righ t -hand side of the 
plot) and the 20 amplitude values from analysing Kilken ny et alj JT999) data (left- 
hand side), mean amplitude values from previous campaigns ar e plotted, appearin g 
as single stars. These are taken from (in chronological order) Koen et al. (1998), 
Ivan Sp aandon k et al.l (|2008), O' TooleetaH J200 0) and Falter et alj (120031) . Note, 
however, that for the sinusoid fit that is also shown in the p l ot the original mean 
amplitude ~ 0.0274 mi for the cam paign by Kilke nny et alj dl999i) was used. In 
the data from |Q' Toole et alj (120001) . f\ could not be found at all. Therefore, the 
amplitude value was set to zero. However, f\ could have had an amplitude bigger 
than zero that was hidden in the noise. This is indicated by the arrow. Since the 
amplitude value sensitively depends on possible unresolved beating or destructive 
interference with noise peaks, values for the amplitude error are not given. 



the main mode, we divided the non-weighted total light curve into 40 chunks, each con- 
sisting of a comparable amount of data points, and fitted the main mode's amplitude and 
phase for each part. The "ghost" frequency f\ a was excluded from the fit, as we consider 
it to be an artefact from imperfectly overlapping data sets (see Sect. 13.31 ) or amplitude 
variation. The result was a linear decreasing trend in the amplitude as well as the phase. 
The latter trend could be fixed by adding 0.00138 cd _1 to the initial frequency guess 
/guess * 179.34699 cd _1 , delivering the corrected value f x * 179.34837 cd~\ while 
the linear amplitude decrease proved to persist at a rate of (-0.3395 ± 0.0089) mmi d _1 . 
In order to put this result into a wider context, we compared the main mode's 
amplitude value and change rate with previous photometric investigations. A plot of 
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the amplitudes versus time is shown in rig. |2] See the caption for the origin of the data 
points. As the simplest reason for amplitude variation is beating between two close 
unresolved frequencies, we tried to fit the data with a sinusoid function. The result with 
a period of ~ 630 d matches the data remarkably well. 



As a further test, we looked for amplitude variation in the data from Kilkenn y et al. 



( 1999), applying the same method as for f\ in the MSST/WET data, but using only 20 
chunks. As fa in those data is closer to f\ than the frequency resolution (~ 0.07 c d" 1 ), 
it was excluded prior to the analysis. Apart from sinusoidal beating between f\ and fa, 
the result was a linear decrease in the amplitude, similar to the MSST/WET data, but 
with a smaller slope. Since the proposed sine indeed expresses a shallower slope there, 
this supports the beating hypothesis. 

Notably, a linear decrease with a slope twice as steep could be found in the 
MSST/WET data for fa, which is the first harmonic of f\. 

5. Multiplet structure 

The fr equency triplet consis ting of fa fa and fa (see tableQ} does also exist in the data 



from iKilkenny et all (1 19991) . however with a different amplitude distribution: Whereas 



the side peaks both had an amplitude of ~ 15mmi before (putting them on ranks "2" 
and "3"), they weakened to ~ 6 mmi in the MSST/WET data. The central peak - now fi 
- was barely detectable as fa with an amplitude of ~ 1 mmi. Furthermore, the splitting 
of the triplet reduced from ~ 0.07 cd" 1 to ~ 0.04 cd _1 , both being in the area of the 
respective frequency resolution. This implies that the true splitting of the triplet is even 
closer and probably has remained constant rather than having changed over five years 
(Suzanna Randall, private communication). 

On the other hand, it is also possible that the triplet struc t ure is an artefact appear- 
ing due to amplitude variation. However. iGuggenberger et alj (12008) demonstrated that 
the triplet (or effectively a doublet due to the weak central frequency) in the data from 



Kilkenny et all ([1999) indeed is real. 



To test if the "new" triplet is real, we excluded f(, and fa from the frequency 
set and checked fi for amplitude variation. The result (see fig. |3j shows a strong in- 
dication for a sine-like periodic function, with both the same period in the amplitude 
and phase. Although beating involving three frequencies is far more complicated than 
beating between two frequencies, the result is clear enough to support the hypothesis. 

In addition to this triplet, other multiplet candidates with a similar spacing could 
be found in the (complete) frequency table, also showing signs indicative of beating 
(see Langfellner et al., in preparation). 



6. Discussion 

6.1. PG 1605+072: a slow rotator? 

The analysis of the triplet comprising f%, fe and fa indicates that these frequencies are 
real, although the exact frequency spacing is not known due to the limited frequency 
resolution. However, at least an upper limit can be given with ~ 0.05 c d" 1 . 

The most plausible explanation for the existence of the multiplets is rotational 
splitting, lifting the m degeneracy. For solid, slow (P rot > 2.5 h) rotat ors and consid- 



ering /j-modes, the frequency spacing is governed by P mt ~ l/Afy ( Charpinet et al. 
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Figure 3. O-C diagram of fi, f(, and /10 when only considering fa. The amplitude 
and phase vary like a sine with the same period, which is indicative of beating. 



2008), including the rotation rate V mt - l/P ro t and the frequency spacing A/#, de- 
pending on the radial order k. Assuming the upper limit for the frequency spacing of 
~ 0.05 cd _1 , a minimum rotation period of about 20 d is inferred. Using crude esti- 
mates and considering both the high-mass and evolved status scenarios, the equatorial 
rotation velocity is then found to be V ea < 0.9 km s -1 or < 1 . 3 km s -1 r espec tively. 

This disagrees with the claims of Heber et al. ( 1999) and Kawaler ( 1999) that have 
been discussed in Sect. CD by orders of magnitude. However, they assumed rotational 



broadening or non-trapped modes. Since Ku assivi et al. 



(2005 ) introduced th e possibil 



O'Toole et ail (l2005h . it is very 



ity of pulsational broadening, which was confirmed by 

possible that PG 1605+072 is a slow rotator. 

Theoretical models by Ivan Spaandonk et all (|2008) and Ivan Grootel et all (1201 Ol) 
cover multiplet spacings of ~ 3.5 cd" 1 or ~ 7.8 cd -1 (~ 40 //Hz or ~ 90/iHz, meaning 
fast rotation) as well as slow-rotation models. Indeed, the reasonable assumption that 
the main mode is radial - as it is part of no multiplet with a spacing around ~ 0.05 c d~ ! 
- constrai ns the available mode l s suffi ciently well to leave a compatible slow-rotation 
model by van Spaando nk et all (l2008h . implying a high mass M ~ 0.7 M Q . Since 



van Grootel et al. (2010) came to the conclusion that all the models they presented pre- 
dicted a mass M ~ 0.76 M Q , this not only supports that PG 1605+072 is a slow rotator, 
but also that it possesses a high mass rather than an evolved status. 
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6.2. Understanding the main mode's amplitude variation 

The sine fitted to the main mode's amplitude data from various campaigns is able 
to describe the positio n of the data points as well as the slope in the data from the 
MS ST/WET and Kilkenny et alj i 19991 ) campaigns, making us believe it to be correct. 
However, the question then arises of why the main mode varies with such a long 
period. For beating, a very close frequency doublet would be required. If the multiplets 
discussed in the previous section are due to rotational splitting, the long-term beating 
must have another reason. If, on the other hand, the long-term beating would be treated 
as being due to rotational splitting itself, extremely low rotation rates in the order of 
0.03 km s _1 would be the result, which seems far-fetched. Should the frequency doublet 
be caused by another physical effect, it is - so far - unclear what specific mechanism is 
applicable. 



7. Summary and conclusions 

The prewhitening of the combined MSST and WET XCov22 photometric data deliv- 
ered 85 identifiable frequencies, with merely small deviations if a weighted light curve 
was used. At least 21 of these frequen cies are either har monics or frequency sums. 
A comparison with the results from Kil kenny et all (Il999h . representing the most ex- 



tensive photometric data on PG 1605+072 prior to MSST/WET, yielded 32 common 
frequencies out of 55. In total, the MSST and WET data increase the number of simul- 
taneously identified frequencies by 30. 

The main mode's amplitude was found to decrease linearly on the timescale of the 
observations. Putting this result in the context of the previous campaigns, we found that 
the amplitude variation can be described by a sinusoid implying long-term beating with 
a period of ~ 630 d. The validity of the sine fi t was supported by co mparing its slope to 
the amplitude variation rate in the data from iKilkenny et all (1 1999b . 



We found that the frequencies with the amplitude ran ks "2", "6" and "10" form 
a real triplet, which has also been present in the data from IKilkenny et all (Il999l) . yet 



changed its amplitude values. Furthermore, the observed frequency spacings have be- 
come smaller, reflecting the better frequency resolution. This indicates that the real 
spacings might be still closer. Additionally, similar multiplet features were found for 
many other frequencies. A more detailed description will be available in Langfellner 
et al. (in preparation). 

The multiplet structure is most likely caused by rotational splitting, indicative of 
PG 1605+072 being a slow rotator. This is in agreement with studies of white dwarfs 
and other single sdB stars. A comparison with theoretical asteroseismic models sup- 
ports the slow-rotation hypothesis and favours the high-mass scenario instead of an 
evolved status. Consequently, this would require an alternative formation channel. 

These results help to further constrain the possible features of PG 1605+072. The 
photometric data and the extensive frequency tables will be published in Langfellner 
et al. (in preparation) and made available online for further investigation. 
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